Very-High-Energy Gamma-Ray Signal from Nuclear Photodisintegration 
as a Probe of Extragalactic Sources of Ultrahigh-Energy Nuclei 



O 



Kohta Murase^'^ and John F. Beacom^'^'^ 
^ Yukawa Institute for Theoretical Physics, Kyoto University, Kyoto, 606-8502, Japan 
^ CCAPP, The Ohio State University, Columbus, Ohio 43210, USA 
^ Department of Physics, The Ohio State University, Columbus, Ohio 43210, USA 
^ Department of Astronomy, The Ohio State University, Columbus, Ohio 43210, USA 

(Dated: July 1, 2010) 

It is crucial to identify the ultrahigh-energy cosmic-ray (UHECR) sources and probe their unknown 
properties. Recent results from the Pierre Auger Observatory favor a heavy nuclear composition for 
the UHECRs. Under the requirement that heavy nuclei survive in these sources, using gamma-ray 
bursts as an example, we predict a diagnostic gamma-ray signal, unique to nuclei — the emission 
of de-excitation gamma rays following photodisintegration. These gamma rays, boosted from MeV 
to TeV-PeV energies, may be detectable by gamma-ray telescopes such as VERITAS, HESS, and 
MAGIC, and especially the next-generation CTA and AGIS. They are a promising messenger to 
identify and study individual UHE nuclei accelerators. 
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I. INTRODUCTION 

The origin of ultrahigh-energy cosmic rays (UHECRs) 
— cosmic rays with energies above 10^*'^ eV — is one 
of the biggest mysteries in astroparticle physics, de- 
spite decades of efforts in theory and experiment [ll, Q ■ 
The Pierre Auger Observatory (PAO) and the High- 
Resolution Fly's Eye (HiRes) are providing newly precise 
clues, though these do not give a consistent picture. 

The UHECR spectrum shows a cutoff around 
60 EeV 0, consistent with the attenuation of 
ultrahigh-energy (UHE) protons due to photomeson 
interactions with the cosmic microwave background 
(CMB), but also consistent with the attenuation of UHE 
nuclei due to photodisintegration interactions with the 
cosmic infrared background (CIB) [5]. In either case, the 
spectrum at the highest energies must be replenished by 
sources within ^ 100 Mpc. The PAO data show a corre- 
lation between UHECR directions and nearby large-scale 
structure @ , which favors protons over nuclei, due to the 
smaller magnetic deflections; the HiRes data do not Q- 

Direct measurements of UHECR composition, through 
the average depth of shower maximum, Xmax, as well as 
the r.m.s. fluctuations around it, are now interestingly 
precise, and thus also probe the identities and properties 
of the sources. This requires comparison to uncertain 
models of hadronic interactions, though the uncertainties 
are much less for the fluctuations. In both quantities, the 
PAO data favor a heavier composition above 10^^'^ eV l8|; 
the HiRes data on Xmax favor a proton composition [9|. 

To reconcile these inconsistent clues, not only UHECR 
observations but also other probes would also be needed, 
e.g., gamma rays and neutrinos that point back to their 
sources. This is especially true if UHECRs are nuclei. 
Their deflections in cosmic or Milky Way magnetic fields 
are larger than than those of protons, so that identifica- 
tion of sources seems more difficuh (e.g., Refs. The 
unprecedented precision of the PAO composition data, 
especially for the fiuctuations around Xmax, motivates us 
to further consider a heavy composition. Though some 



heavy composition models have been suggested (e.g., 
Refs. 0), most papers about gamma-ray and neutrino 
signals have focused on proton sources. 

Very-high-energy (VHE) gamma rays and neutrinos 
can be produced by sources of UHECR nuclei as well 
as protons. However, nuclei are special because they can 
produce nuclear de-excitation gamma rays following pho- 
todistintegration interactions. We show that this pro- 
vides a characteristic energy spectrum that could differ- 
entiate extragalactic accelerators of nuclei and protons. 



II. VHE GAMMA-RAY AND NEUTRINO 
PRODUCTION 

The most widely-discussed UHECR sources are 
gamma-ray bursts (GRBs) (Til [l^ and active galac- 
tic nuclei (AGN) [li [11. The predicted VHE spec- 
tra of gamma rays and neutrinos can be calculated with 
some assumptions, e.g., the yields from photomeson in- 
teractions depend on the cosmic-ray spectrum and tar- 
get photon spectrum, and can strongly vary between 
source models (see, e.g., Refs. [isl - lisl for GRBs and 
Refs. [L9-21] for AGN). GRBs and AGN can also be con- 
sidered as sources of UHE nuclei [H, [H, [H, [l^ . In this 
work, as a demonstrative example, we mainly consider 
GRBs [ll, and especially low-luminousity (LL) bursts 
such as GRB 060218, which might be more numerous 
than classical GRBs [25|]. Applications to other sources 
such as AGN would be possible for given source models. 



We assume a cosmic-ray energy spectrum of 
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where E"^""^ 



is the maximum energy. We 
adopt p — 2.3, consistent with observations iS]. The 
normalization is determined from the UHECR energy in- 
put rate, estimated as E"^^^^^ ~ 10'''* erg Mpc'^yr"* at 



Eq — 10*® eV (e.g., Refs. jl|). For persistent sources, this 
is the product of the UHECR luminosity per source and 
source density. For transient sources, this is the product 

of the UHECR energy per burst 



£;2£[^ and the 



dE 



2 



apparent rate p. GRBs with p ^ 1 Gpc '^yr ^ require 
E^'^CK = ^crIbo ~ 10^^ erg LL GRBs, hypernovae 
and AGN flares, which may have p~ 10^^"^ Gpc^'^yr^^, 
require E^lcK ^ erg 

For the target photon spectrum in the source, we use a 
broken power law, dn/de cx e""'-'', as for the synchrotron 
emission mechanism. Here e is the photon energy in the 
comoving frame, £ob(~ Te) is that in the observer frame, 
and r is the bulk Lorentz factor. For GRB prompt emis- 
sion, ai_21 1 for e < and at ^ for < e are typical 
values [2J|, and is the break energy. 

Both protons and nuclei can undergo photomeson in- 
teractions above the pion production threshold. For nu- 
clei, the cross section and inelasticity are ~ Aup^ 
(neglecting shadowing) and n^^s ~ i^p-^/A [12, 1201 ■ Using 
the A-resonance approximation, the photomeson produc- 
tion efficiency /mes ~ tdyn/tmcs is (e.g., 
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'rl is 

the resonance energy, L?, is the photon luminosity at e^^, 



where E 
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1.8 X 10^^ eV {A/56){e 
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is the emission radius, td 



yn w r/Fc is the dynamical 
time scale of the relativistic source, and tmcs is the pho- 
tomeson energy loss time scale. Here, multi-pion produc- 
tion, which is important at sufficiently high energies, is 
not considered (c.f. Eq. (18) in Ref. [12]). 

The value of /mes has large uncertainties coming from 
parameters such as r. As in previous works [H, [Tsl - 
1201, one can consider cases where UHE nuclei are effi- 
ciently disintegrated in the source, i.e., where /mcs ~ 
0.01 — 0.1. But in this work, since we consider the PAO 
case where UHECRs are largely heavy nuclei, we focus 
on the scenario where GRBs are sources of UHE nuclei. 
As shown in Refs. [l^, H^Ij besides several open issues 
such as particle escape, UHE nuclei rather than UHE 
protons are produced and survive in both the prompt 
and afterglow phases, when r is large enough. We con- 
sider such cases where UHE iron can completely sur- 
vive photodisintegration {ta^ ^ 1) [H, 113, 1131, where 
/mcs S 1-9 ^ 10~^(A/56)~^'^^, and note that our purpose 
here is not to propose GRBs as sources of UHE nuclei. 

Nuclei interact with target photons via the photodis- 
integration interaction, A + j ^ A'* + X [2a]- The 
main contribution comes from the giant dipole resonance 
(GDR), which decays by the statistical emission of a 
single nuclcon, though at higher energies, quasi-deutron 
emission and fragmentation occur. After photodisinte- 
gration, the daughter nucleus is typically left in an ex- 
cited state, which can immediately emit gamma rays as 
A'* A' + j. The gamma-ray multiplicity is ~ 1 — 3 and 
the energy in the nuclear rest frame is ~ 1 — 5 MeV. We 
adopt Epe = 2.5 MeV and npe = 2 [13], so that the energy 
fraction carried by gamma rays is Kdeex ~ 10~^(56/A). 

The de-excitation process is not often taken into ac- 
count for gamma-ray production. Examples of applica- 



tion to Milky Way sources include Refs. [28|,l29|. Here we 
consider this process in an extragalactic context, which 
is relevant because we require that UHE nuclei from 
nearby sources survive sufficiently to reproduce the heavy 
UHECR composition deduced from the PAO data. 
In the GDR approximation, the dc-excitation efficiency 

^/rieex « tdyn/tdccx ~ Kdcoxtdyn/tyl7 beCOmCS (e.g., [22, 
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is the resonance energy and tA-^ is the photodisinte- 
gration time. The observed energy of de-excitation 
gamma rays from iron will be e'^'^'^^^ « T^a^Fc — 
1.5 TeV {A/b&y'^EA,i6.i,. Note that we have 

/mcs(4"r^) - 20(A/56)-°'Vdccx(41^^)- We can see 
that the de-excitation efficiency is typically low when the 
nucleus-survival condition is satisfied, and the photome- 
son process is more efficient. Nevertheless, we see that 
the de-exicitation process leads to the interesting signal 
for nearby sources, and its component may dominate over 
the other components around the TeV scale. 

Gamma rays may be attenuated in the source by 
pair creation. However, when the conditions are such 
that UHE iron survives in the sources (i.e., rFc7 < 1), 
roughly speaking, we can neglect gamma-ray attenua- 
tion. This is because a-po-y — 8 x 10~^^ cm^ O.ltiT ~ 
a~f^ so that t-^^ ^ Tpc-y- (Of course, note that the de- 
tail depends on the target photon field 12]). Gamma- 
ray attenuation also happens en route, even for the close 
sources we consider, and this is discussed below. 

Neutrinos are also produced by nuclei. The typical 
energy of conventional photomeson neutrinos is eI""^^"^ ~ 
0.28 PeV {A/56y Ea^7.5- Neutrinos also come from the 
/3-decay of nucleons from photodisintegration. A nucleon 
carries kgdr ~ 1/A of the energy of the nucleus, and 
about half of nucleons are neutrons. The neutrino energy 
in the neutron rest frame is ~ 0.48 MeV [1^, so that 

~ 0.29 TeV {A/56)~^Ea- 
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III. 
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PREDICTED GAMMA-RAY AND 
NEUTRINO FLUENCES 



In the previous section, we analytically evaluated the 
nuclear de-excitation efficiency, and we can predict spec- 
tra of de-exicitation gamma rays from Eq. (2). In this 
work, however, we perform much more elaborate calcu- 
lations, using detailed cross sections of photodisintegra- 
tion, photomeson, and Bethe-Heitler processes [16,, ^3(1 , 
where various cooling processes are taken into account. 
Following the method used in Refs. [ll,lil], the photodis- 
integration time and other cooling time scales for nuclei 
and mesons (photodisintegration t^is, photomeson imcs, 
Bethe-Heitler teH, synchrotron i^syn, inverse-Compton 



3 



- 


1 1 


^•^ 1 1 


1 




t-L t-' 

Lad , Ldyn 


Lace -.^ 




- 




Ld^ . — 








■' ' ' 

''" tAsyn 


^deex 
1 




2 4 


6 8 


10 



log(E^/r [GeV]) 

FIG. 1: The various cooling time scales and acceleration time 
for iron in the internal shock model for LL GRBs. Time 
scales are measured in the comoving frame of the outflow. 
The source parameters are r — 10^^'* cm, F — 10, l}^ — 
lO^'*'^ erg s-\ a.1 = 1 and an = 2.2, and e^t = 5 keV [3. 
Note that the disintegration loss is usually dominated by emis- 
sion of nucleons, i.e., t^l^ is larger than t~^^^^^ tZ\as and tg^. 
Here, fdis is the photodisintegration energy loss time, which is 
longer than the photodisintegration interaction time Ia-i, and 
tdccx is the de-excitation energy loss time due to gamma-ray 
emission. Note that tdocx ~ ^bh at high energies, suggest- 
ing the potential importance of the de-excitation process. In 
fact, the de-exicitation process leads to production of higher- 
energy gamma rays, so that its component is more important 
at the VHE range in our case (see text). 



tic, E^nd adiabatic loss tad) are numerically evaluated 
(see Fig. 1), and emissions by photomeson and cosmic- 
ray synchrotron processes are also calculated in detail. 
The magnetic field is determined from the equipartition 
between the magnetic energy density and the target pho- 
ton energy density, which is often assumed in the classical 
optically synchrotron scenario for GRB prompt emission. 

The maximum energy E^^^ is estimated by comparing 
the acceleration time tacc ~ '^'^'z^BcT ^^^^ with the dy- 
namical and cooling times, and E'^'^'^ ~ {Z/ 26)10'^° eV 
is obtained for our parameters described in the caption 
of Fig. 1. (Note that tacc = zeBcr ^^^^ E3- 
We use a more conservative expression in this work.) 

We also use approximate formulae in Refs. 29] to eval- 
uate primary spectra of gamma rays and neutrinos from 
photodisintegration. In addition, the gamma-ray attenu- 
ation by pair creation with the CIB/CMB is considered. 
The mean free path of gamma rays for pair creation and 
energy loss length of electron-positron pairs for inverse- 
Compton are shown in Fig. 2. 

Motivated by the PAO results, we consider a heavy 
composition of 50% iron and 50% proton (to be conser- 
vative, not 100% iron). The results for the de-excitation 
gamma rays would be unchanged as long as the composi- 
tion is largely heavy nuclei of some species. The condition 
TF07 < 1 and < 1 is satisfied for our adopted source 
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FIG. 2: The mean free path of high-energy photons for pair- 
creation and the energy loss length of electron-positron pairs 
for inverse- Compton in the Universe, respectively (solid and 
dashed lines). The dotted curves show the synchrotron cool- 
ing length for given IGMF strengths. 



parameters [Hj]. Also, the corresponding baryon load- 
ing factors are 

^cr/4'° ^ 10^"^. Although such large 
values might be possible, smaller values which might be 
preferred can also be achieved for flatter or broken power- 
law cosmic-ray spectra [l^, [13] • The origin of the nuclei 
and the baryon loading depend on source physics and 
environments (e.g., Refs. ]12, 22] for GRBs). 

The resulting gamma-ray fluxes are shown in Fig. 3. 
At the highest energies, the primary photomeson gamma 
rays are dominant; however, these are severely atten- 
uated by the CIB/CMB. Detections of the cosmic- 
ray synchrotron component (peaking at Ey'^^^'^^ ~ 
33 MeV (Z/26)(A/56)-3£:2 ^gBi gFJ-i) is also diffi- 
cult, due to its insufficient fiux. (The proton compo- 
nent typically seems more important in cosmic-ray syn- 
chrotron emission at higher energies, since e!^'^^^'^^ oc 
E\ZA~^ [13)- In any case, not only nuclei but also pro- 
tons, contribute to those meson and electromagnetic pro- 
duction processes. Hence, even if signals were detected, 
it would be difficult to prove that nuclei are accelerated 
in extragalactic sources. 

The de-excitation process is invaluable as a probe of 
UHE nuclei accelerators. As in Eq. (2), it has a charac- 
teristic spectrum in the VHE range, reflecting the tar- 
get photon spectra, though the details depend on source 
models. This allows it to be distinguished from other pro- 
cesses, once the target photon spectra are known. The 
VHE range is doubly favorable: the attenuation due to 
the CIB/CMB is modest for nearby sources, and this is 
where the sensitivity of gamma-ray telescopes is best. 

The fluence of de-excitation gamma rays shown 
in Fig. 3 can be reproduced with a simple es- 
timate. Using Eq. (2), we obtain £'^(^7'^°'^'^'' 
(/docx4R,Fo)l£<Fo.)/(47ri:)2) _ 5 X 10-8 erg cm-^ at the 
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FIG. 3: Energy fluences of gamma rays from a LL GRB with 
^HECR = 10'^° '^ erg and distance D = 100 Mpc. The source 
parameters are indicated in the caption of Fig. 1 with emission 
duration T — 5000 s. The red curves are the de-excitation 
component with attenuation (thick solid line), without atten- 
uation (thin dotted line), and with a possible cascade com- 
ponent (thin dashed line). The blue curves are the same 
for the photomeson component, and the green curve is the 
ion-synchrotron component. For comparison, assuming that 
prompt emission is synchrotron, the SSC component is also 
shown (yellow curve). The sensitivities of Fermi/LAT 'sT] 
and VERITAS (with a duty factor of 20 %) [32] are labeled. 



r^. Although we 



break E^r""^ 0.5 TeV (£^^,.5 kcv) 
have used typical parameters inferred from the observa- 
tion of a LL GRB, one can explicitly see the parameter- 
dependence of the de-excitation gamma-ray flux from 
Eq. (2). For GRBs at D < 100 Mpc, the de-excitation 
signals might be detected by present-generation gamma- 
ray telescopes such as VERITAS (Very Energetic Ra- 
diation Imaging Telescope Array System), HESS (High 
Energy Stereoscopic System), and MAGIC (Major At- 
mospheric Gamma-ray Imaging Cherenkov Telescope). 
The prospects with next-generation telescopes such as 
CTA (Cherenkov Telescope Array) and AGIS (Advanced 
Gamma-ray Imaging System) are even more encourag- 
ing. HAWC (High Altitude Water Cherenkov observa- 
tory) (with TeV sensitivity of 5 x 10~^ erg cm~^ for 
T — 5000 s) would also be useful because of its very large 
field of view and duty cycle. 

When gamma rays are attenuated by the CIB/CMB, 
the secondary electron-positron pairs created will up- 
scatter CIB/CMB photons to gamma-ray energies, and 
the process repeats until the pair creation threshold is 
reached. The detectability of this secondary cascade 
emission depends critically on the intergalactic magnetic 
field (IGMF) . Especially, since the mean free path of TeV 
photons is hundreds of Mpc, pairs with ~ TeV energies 
are likely to feel the IGMF in voids, which is very un- 
certain. (On the other hand, UHECRs should also feel 
the stronger IGMF in the structured region (clusters, fil- 
aments and sheets) and Milky Way magnetic field [3J].) 



FIG. 4: Energy fluences of neutrinos from a LL GRB, as in 
Fig. 3. For neutrinos from neutron decay, the time-integrated 
(thick solid line) and early (within T = 5000 s; thin solid line) 
cases are shown. Neutrinos from pion decay are also shown. 
Neutrino mixing is not taken into account. 



If the IGMF in voids is strong enough (Big > 10"^*^ G), 
then the emission becomes nearly isotropic and its flux 
is diminished due to the magnetic time spread, so the 
cascaded gamma rays are not detectable [3a |. 

However, if the IGMF in voids is weak, there may 
be a detectable pair echo for transients or a pair halo 
for persistent sources [H, Is^]. The echo duration is 
AtB ~ 10'* s E~^2BiG -20' where E^, is the energy of 
pairs The secondary cascaded gamma-ray spec- 

tra are shown in Fig. 3, which are calculated followin, 
Ref. [l^. (Note that the calculations shown in Ref. [l 
were properly performed only for energies above the TeV 
scale, due to the focus on UHE gamma rays. Hence, the 
interaction and attenuation lengths of gamma rays and 
pairs in the Universe, shown in Fig. 1 of Ref. [18], are cor- 
rect above the TeV scale.) Thus, for a sufficiently weak 
IGMF, we expect that the cascaded photomeson gamma- 
ray flux, overwhelming the de-excitation gamma-ray flux, 
is detectable. However, it may be possible to discrimi- 
nate between non-cascade and cascade components using 
temporal information 35|. 

Depending on source models, there are possibilities for 
which the de-excitation component is not dominant in 
the VHE range. There would be synchrotron from pairs 
from decaying muons (originating from the photomeson 
process) and pairs from the Bethe-Heitler process. 
In GRB cases, photomeson production is typically 
more important than the Bethe-Heitler process at high 
energies (/mcs > /bh ~ 10~*), and its energy fluence 

' " /8/{4t:D^) with 
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peak of 4^''"^ 
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The estimated fluence at TeV is '--^ 2 x lO^'* 
for the pure iron case. Noticing 
7.6 X 10-31 cm2(ZV^) at e 



erg cm 

'^bhO'bh ~ 
10(2TOec2) [30], the 
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energy fluence of synchrotron emission by pairs 
from the Bethe-Heitler process is similarly esti- 

mated as ^ (/bh^DI^cbh, (£;^/4r"^)^^/(47ri?')> 
where the effective optical depth for the 
Bethe-Heitler process is /bh — 0.76 x 

and 4f ) c 5.4 x 10^' eV (A/56)(e^,^5 kcv)"'r2. 
The synchrotron peak is at E^^^'^'^'^ ^ 1.3 x 

10^ eV {A/56y^{Ef^^^^j^fBi,5r^\ and the esti- 
mated fluence at TeV is ^ 6.3 x 10"^ erg cm^-^ for 
the pure iron case. The inverse- Compton emission by 
pairs, which is suppressed at high energies due to the 
Klein- Nishina (KN) effect, is also not dominant at the 
TeV scale. In our GRB cases, a; ~ 1 and p ^ 2.3 lead to 
detectable de-excitation signal, and we may expect that 
it dominates over the photomeson and Bethe-Heitler 
signals at the TeV scale. Generally speaking, however, 
the relative importance of the de-excitation signals 
depends on B, a, p, and the composition. 

There may also be a purely leptonic component such 
as synchrotron self-Compton (SSC) emission. However, 
in our cases, the KN effect leads to significant suppres- 
sion of SSC emission at very high energies. If keV-MeV 
emission is synchrotron from fast-cooling electrons, from 
e^jj ~ keV and Ub = U^, the typical Lorentz factor of 
electrons is ^^,1 ^ 10'*- The KN effect becomes important 
at £;fKN « r7e,.mec2 ^ 50 GeV (< ^ « ^^'eAJ. 
Using the SSC- Y parameter in the Thomson limit, 
y„oKN ~ {-l + ^l + 4{Ue/UB))/2 ~ 0.6, the SSC energy 
fluence at = 1 TeV is estimated as [13], E^c/y^j'^^ « 

2 X 10^* erg cm^^ (see Fig. 3). Note that similar sit- 
uations where the KN effect is significant are expected 
for GRB 080916C [S^. But this is expected when the 
typical Lorentz factor of electrons is large enough, and 
resulting fluence depends on parameters such as B, T, 
'^HECR ^^'^ unknown GRB radiation mechanisms (which 
affect a). For example, if je,i ^ 10"*, F ^ 500, and 
EqJj ^ 250 keV (which are expected for high-luminosity 
GRBs), the KN suppression occurs at 2.5 TeV, while 
the de-excitation component peaks at ^ 25 TeV and 
may be visible around this VHE range. Therefore, 
the de-excitation and the other hadronic signals can 
dominate over the purely leptonic components at the 
VHE range, when the KN effect is significant. In our 
GRB cases, we expect that the de-excitation component 
may be visible above the TeV range (see Fig. 3), but the 
de-excitation and pionic gamma-ray components seem 
to be overwhelmed by electromagnetic (electron- and 
ion-synchrotron) components around the GeV range. 

Then, at least in our cases, the "hard" de-excitation 
component is visible at > TeV over those from other pro- 
cesses in the sources or en route, though variations within 
the model uncertainties might change this. The purpose 



of this work is to demonstrate the potential importance 
of the de-excitation signals from UHECR sources, and 
to motivate VHE gamma-ray searches. Hence, here we 
avoid further studies, and expect that uncertainties will 
be reduced once sources are measured by future searches. 

^In Fig. 4, the neutrino fluences are shown. For Ice- 
Cube detections, fluences of > 10~^eI^^^ erg cm"^ at > 
0.1 PeV are required [lH|. The time coincidence is impor- 
tant to reduce atmospheric neutrino backgrounds, but it 
would not work well for neutrinos from photodisintegra- 
tion. Relativistic neutrons have lifetimes of i„ ~ ja 887 s 
in the comoving frame, so that the fluence is suppressed 
at energies above £;^°'" ~ 0.27 GeV (>l/56)rfT3.7. Thus, 
when nuclei survive efficiently, we expect that their 
source neutrino fluxes are too low to be detected. 



IV. SUMMARY AND DISCUSSIONS 

We calculate the gamma-ray emission from UHE nu- 
clei sources, inspired by the surprising FAQ data sug- 
gesting that the UHECRs are largely heavy nuclei. We 
show that the de-excitation process of gamma-ray pro- 
duction, which is unique to nuclei, could be distinguished 
from other gamma-ray production processes and could 
be detected by VHE gamma-ray telescopes. Therefore, 
it could identify extragalactic accelerators, with the de- 
tectable gamma rays probing cosmic-ray nuclei at the low 
end of the UHE range. Importantly, our results show 
the most promising technique for identifying UHE nu- 
clei sources, as the cosmic-ray nuclei will be strongly 
deflected [lO|, the neutrino fluxes are low [l^, HI], and 
other hadronic gamma-ray signals such as pionic gamma 
rays cannot be cleanly identified. If the nuclear sig- 
nals are seen from UHECR sources, that would favor the 
PAO claim of a heavy nuclear composition. If they are 
not seen, then further measurements and studies will be 
needed to know the nature of UHECR sources. Obvi- 
ously, all of the cosmic-ray, neutrino and gamma-ray ex- 
periments will be important (e.g., detecting UHE gamma 
rays [3 and/or cosmogenic neutrinos [39j] would favor 
protons over nuclei). 

Detections would be limited to nearby and/or energetic 
sources, because (1) the de-excitation efficiency is very 
low when iron nuclei survive, and (2) TeV-PeV photons 
from distant sources are attenuated by the CIB/CMB. 
If UHECR sources are transient, the event rate within 

100 Mpc is - 1.3 (^IIecr.so.s)"^ y^"^ -18]. These events 
can be found by low-energy all-sky monitors and followed 
up with VHE gamma-ray telescopes. 

Visibility of the de-excitation component may also 
be limited by other competing processes. As discussed 
in this work, the de-excitation process is important, 
but other photomeson, Bethe-Heitler, and purely lep- 
tonic processes such as SSC emission may be considered. 
Purely leptonic emission would be likely to be most im- 
portant. If the baryon loading is not small and the KN 
suppression is significant, the hadronic signature may be 
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observed at very high energies. Very high-energy pairs 
from the photomeson and Bethe-Heitler processes also 
lead to synchrotron and inverse-Compton emission. For 
a ~ 1, the Bethe-Heitler energy loss rate of nuclei is 
roughly comparable to the de-exiction energy loss rate, 
but the typical energy of the de-excitation gamma rays 
is higher than that of gamma rays from pairs, which 
may make the de-excitaton signal observable (see above). 
Note that the photomeson energy loss rate is larger than 
that of the Bethe-Heitler energy loss rate at high energies 
in our cases (see Fig. 1), but the Bethe-Heitler process 
can be more important if photon indices are steeper (e.g., 
a - 2). 

Despite these limitations, the most attractive feature 
of the nuclear de-excitation signal is that it is unique 
to nuclei and cannot be produced by protons. On the 
other hand, the photomeson and Bethe-Heitler processes 
are also induced by protons. The nuclear Bethe-Heitler 
process leads to more gamma rays than the proton Bethc- 
Hcitler process for a spectrum with given normal- 
ization. However, discrimination between protons and 
nuclei is not easy, unless we know each normalization 
of spectra of proton and nuclei. Generally speaking, all 
those three processes (de-excitation, photomeson, and 
Bethe-Heitler processes) can have important roles in 
UHE nuclei sources. 

Our calculations on the de-excitation signal are poten- 
tially important and formulae (see Eq. (2)) are general. 
We have performed the detailed calculations for prompt 
emission from GRBs, but detectability is different for 
other models and associated parameters. For example, 
one can apply our formulae to the external shock model 
for afterglows. Although GRBs have been studied in this 



work as a detectable example, it would be possible to con- 
sider this process for other candidate UHECR sources, 
such as AGN and hypernovae. As for AGN, though vari- 
ous acceleration and emission zones (e.g., blazar regions, 
hot spots and cocoon shocks) can be considered, the con- 
dition of nuclei surviving requires large emission radii, 
as discussed in Ref. [2^. This means low photomeson, 
Bethe-Heitler, and de-excitation efficiencies, which limit 
our accessibility to UHE nuclei sources. In fact, the nu- 
cleus survival condition taj < 1 @| gives /mcs ^ 1-9 x 
10"3(A/56)"^-^\ /bh < 2.6 X 10-4(Z/26)^(A/56)"^-^^ 
and /docx < 0.95 x 10"^(A/56)"\ In the case of blazars, 
we typically expect strong SSC and/or external inverse- 
Compton emission, which may mask those weak hadronic 
signals. Radio galaxies with radio lobes and /or hot spots 
may also be viable candidates of nuclei [IJ]. Low ef- 
ficiency of gamma-ray and neutrino production implies 
that detections of the nuclear signals may be possi- 
ble only for nearby sources. In this sense. Gen A (at 
^ 3.8 Mpc) may be of special interest. Some authors 
argued that all the UHECRs may come from Cen A as 
heavy nuclei [40| . Such the single source scenario can 
potentially be tested by constraining existence of the de- 
excitation and the other hadronic components. 
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